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Abstract 

We present the results of a study of heavy-light-quark bound states in the 
context of the reduced Bethe-Salpeter equation with relativistic vector and 
scalar interactions. We find that satisfactory fits may also be obtained when 
the retarded effect of the quark-antiquark interaction is concerned. 
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Because of the limited understanding for confinement at present, more theoretical efforts 
to be made related to this issue are worthwhile. In an earlier paper JTJ, we presented results 
of a relativistic analysis of the spectrum of light- and heavy-quark-antiquark bound states 
based on the reduced Bethe-Salpeter(BS) equation, while retardation effects in the quark- 
antiquark interaction kernel were taking into consideration. The results are stimulating and 
appears having clarified the problem pointed out by Durand et al. [0 for the static scalar 
confinement in reduced Salpeter equation. The "intrinsic flaw" of the Salpeter equation with 
static scalar confinement could be remedied to some extent by taking the retardation effect 
into confinement. In the on-shell approximation for the retardation term of linear confine- 
ment, the notorious trend of narrow level spacings for quarkonium states especially for light 
quarkonium states is found to be removed. A good fit for mass spectrum of S-wave heavy 
and light quarkonium states (except for the light pseudoscalar mesons) is obtained using 
one-gluon exchange potential and the scalar linear confinement potential with retardation 
taken into account. 

In this paper we extend our previous study to the heavy- light- quark systems (Qq or qQ) 
in order to get a fully understanding of the retardation effects on meson spectra. Based 
on the same procedure taken in Ref. [IJ, by solving the reduced Bathe-Salpeter equation 
numerically, the previous conclusion is further substantiated through this study. 

We assume the confinement kernel in momentum space taking the form 

1 _ 1 
-P 2 ) 2 (p 2 ~ Po) 



as suggested by some authors as the dressed gluon propagator to implement quark con- 
finement 0. Here p is the 4-momentum exchanged between the quark and antiquark in a 
meson. If the system is not highly relativistic we may make the approximation 



G(P)«7S— ^«75^ 1 + ^ . (2) 




and may further express po in terms of its on-shell values which are obtained by assuming 
that quarks are on their mass shells. This should be a good approximation for cc and bb 
systems. However, to get a qualitative feeling about the retardation effect considered here, 
we will also use (fj) for heavy-light-quark mesons, though the approximations are not as 
good as for heavy-heavy-quark mesons. With the above procedure, the scalar confinement 
kernel becomes instantaneous again but with some retardation effects have been taken into 
the kernel. In the static limit, the retardation term vanishes and the kernel returns to 
G(q) oc 7^5, which is just the Fourier transformation of the linear confining potential. In 
this paper, we will use this modified scalar confining potential in which the retardation 
effect is incorporated and the one-gluon-exchange potential in the framework of the reduced 
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Salpeter equation to study the mass spectrum of qQ mesons, and the structure of the 
hyperfine splittings of the heavy-light-quark systems will also be investigated in this paper. 

In quantum field theory, one of the basic descriptions for the bound states is the Bethe- 
Salpeter equation ||. We can define the BS wave function of the bound state | P) of a quark 
ip(xi) and an antiquark ^(x 2 ) as 

X{xi,x 2 ) = (0 | Ti/)(x 1 )^(x 2 ) | P). (3) 

Here T represents the time-order product, and the wave function can be transformed into 
the momentum space, 

Xp(q) = e~ lP - x J d'xe-^xixi, x 2 ), (4) 

where P is the four-momentum of the meson and q is the relative momentum between quark 
and antiquark. As the standard measure in solving the BS equation, we choose center-mass 
and relative coordinates as variables, 

X = TJiXi + 772X2, x = xi-x 2l (5) 

where r\i = m ^! m2 (i = 1,2). After taking the Fourier transformation, in the momentum 
space the BS equation reads as 

( J f l -m 1 ) Xp ( q )( J f 2+ m 2 ) = ±- fd 4 kG{P, q -k) X p{k), (6) 

Z7T J 

where p% and p 2 represent the momenta of quark and antiquark, respectively, 

pi = 771 P + q, p 2 = r] 2 P - q, (7) 

and G(P, q — k) is the interaction kernel which acts on x and is determined by the interquark 
dynamics. Note that in Eq.(|6]) vii\ and m 2 represent the effective constituent quark masses so 
that we could use the effective free propagators of quarks instead of the full propagators. This 
is an important approximation and simplification for light quarks. Furthermore, because of 
the lack of a fundamental description for the nonperturbative QCD dynamics, we have to 
make some approximations for the interaction kernel of quarks. In solving the Eq.(^|), the 
kernel is taken to be instantaneous, but with some retardation effect being taken into it; the 
negative energy projectors in the quark propagators are neglected, because in general the 
negative energy projectors only contribute in higher orders due to M—E\—E 2 -C M+Ei+E 2 , 
where M, Ei, and E 2 are the meson mass, the quark kinetic energy, and the antiquark kinetic 
energy respectively. Based on the above assumptions the BS equation can be reduced to a 
three-dimensional equation, i.e. the reduced Salpeter equation, 
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( P o _ Ei _ E 2 )<$> p (q) = A^ 7 °| d 3 kG(P, q, k)§p(k)j°A 2 _. (8) 

Here 

$-(g) = Jdq XP (q ,q), (9) 
is the three dimensional BS wave function, and 

A+ = ^(Ei + 7°7 " Pi + m l7 °), (10) 

A 2 _ = -L(P 2 - 7 ° 7 ■ p 2 - m 27 °), (11) 

are the remaining positive energy projectors of the quark and antiquark respectively, with 



Ei = \jm\ + p{ , E 2 = \jm\ + p 2 . The formal products of G$ in Eq.(^) take the form 

G$ = ]T G.O^O, = G s $ + 7m ® l^Gv^i (12) 

j 

where = 7a4 corresponding to the perturbative one-gluon-exchange interaction and 0, = 1 
to the scalar confinement potential. 

From Eq.(^[) it is easy to see that 

Al®p(q) = $p(q), $p(q)M = $p(q). (13) 

Considering the constraint of Eq.(|TJ) and the requirements of space reflection of bound 
states, in the meson rest frame (P = 0) the wave function $p(q) for the 0~ and 1~ mesons 
can be expressed as 

$°T(g) = A^ 7 °(l + 7 ) 757 °A^(g), (14) 

SjT (?) = A^ 7 °(l + 7 °) ^ 7 °A 2 /(g), (15) 

where ^ = 'y lJj e fl is the polarization vector of 1~ meson; <p{q) and f(q) are scalar functions of 
q 2 . It is easy to show that Eqs.(|TJ]) and ( |i~5| ) are the most general forms of the S-wave wave 
functions for 0~ and 1~ mesons in the rest frame. 

Substituting Eqs.(0), (0), and ([15D into Eq.®, one can derives the equations out for 
<p(Q) and /(?) in the meson rest frame [0, 
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M(fx(q) = (E 1 +E 2 )< fl {q) 

^2 

"d 3 k(G s (q, k) -4G v (g, k))<pi(k) 



-±2 

E X E 2 + 771x1712+ Q 



AExE 2 



AE\E 2 J E\ + E 2 

Ei + E 2 f 3 ^ ^ N mi+m 2 



d d kG s (q,k)(Q • fc) - 1 ' z yi(fc) 



J 'E x m 2 + E 2 m x 

mi-m 2 f 3 - - - - E x -E 2 



with 



/^n (mx +m 2 + + J E 2 )(^i^2 + ^2^1) 
*" ( '" = 4^( m , +ma) * (8) " (17) 



and 



Mf 1 (q) = {E 1 + E 2 )f 1 {q) 

~^e\e 2 I d ^ Gs ^^ *) " 2GV ^' *0)(£i™2 + SamO/iU) 

^1 + ^2 /", 3 ,^ ^^E l 77l 2 + E 2 77l l 



d d kG s (q,k) -f^—FT fi(k) 



AE\E 2 J E\ + £/ 2 

j-2 



<FMG S (?,fc)+4GV(?,fc))(? • fc)/i(jfc) 



— mim 2 + 9 

- ^"^T 1 / *) - 2Gy(q, k))d ■ k)^^hCk) 

E 1 + E 2 -m 2 -m 1 f 3 - - - -\ 2 1 



d 6 kG s {q,k){q ■ k) 2 -———, ; h{k) 



2E l E 2 q 2 J -/ El + E2 + mi + m2 - 

/ d 3 *GV(?, *)(? ■ fc) 2 ■ * ■ „, (18) 



with 



/,(?) = - mi+ 7^ 1 + £2 /(?). (19) 

Eqs.(|i~6l) and (|i~8f ) can also be formally expressed as more compact forms 

[M — Ex — E 2 ) Vl (q) = [d 3 kJ2 k)G { (q, k)^(k), (20) 

J i=S,V 

(M — Ex — E 2 )h(q) = fd 3 k £ ^"(g, k)Gi(q, k)h(k). (21) 
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In case of taking the nonrelativistic limit for both quark and antiquark, expanding in terms 
of Q /mi and Q /m 2 , it can be approved that Eqs.(|T^) and flT8| ) are identical with the 
Schrodinger equation to the zeroth order, and with the Breit equation to the first order. 

To solve Eq.(|D one must have a good command of the potential between two quarks. 
At present, the reliable information about the potential only comes from the lattice QCD 
result, which shows that the potential for a heavy quark-antiquark pair QQ in the static 
limit is well described by a long-ranged linear confining potential ( Lorentz scalar Vs ) and 
a short-ranged one gluon exchange potential ( Lorentz vector Vy ), i.e. ||[7| 

V(r) = V s (r)+ lfl ®YV v (r), (22) 



with 



Vo(r) 



Xr- 



ar 



(23) 



Vy(r) = -tS^e-ar (24) 
6 r 

where the introduction of the factor e~ ar is not only for the sake of avoiding the infrared(IR) 
divergence but also incorporating the color screening effects of the dynamical light quark 
pairs on the "quenched" QQ potential ||. Although the lattice QCD result for the QQ po- 
tential is supported by the heavy quarkonium spectroscopy including both spin-independent 
and spin-dependent effects [0-11]. we will employ this static potential below to the heavy- 



light-quark systems as an assumption. The interaction potentials in Eq.(|22|) can be trans- 
formed straightforwardly into the momentum space, where the strong coupling constant 

, \ 12tt 1 
a.{P) = — 3—. (25) 

QCD 

— ^2 

is assumed to be a constant of 0(1) as V — ^ 0. The constants A, a, a, and Aqcd are the 
parameters that characterize the potential. 

In taking the retardation effect of scalar confinement into consideration, as discussed in 
Ref. Q, the confinement will be approximately introduced by adding a retardation term 
^ to the instantaneous part j-^y as given in Eq.(Q), and p% will be treated to take its 
on-shell values which are obtained by assuming that the quarks are on their mass shells, 
which means that the retardation term will become instantaneous rather than convoluted. 
By this procedure the modified scalar confinement potential will include retardation effect 
and become 
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a 7T 2 (p 2 + a 2 ) 2 



+ 2 ^ (P + k) 2 + m 2 -^k 2 + m 2 ) 2 , (26) 

tx a (p z + a A y 

which is related not only to the interquark momentum exchange p but also the quark mo- 
mentum k itself. 

Based on the formalism and discussions above, we can now embark on the numerical 
calculations, in which we take the following values for input parameters 

A = 0.21 Gel/ 2 , a = 0.06 GeV, a = e = 2.7183, A QCD = 0.19 GeV, C = -0.05 (27) 
and 

m u = md = 0.35 GeV, m s = 0.5 GeV, m c = 1.68 GeV, m b = 4.925 GeV, (28) 

which fall in the scopes of customarily usage. The numerical results with retardation are 
listed in Table I. For the convenience of comparisons, results obtained without retardation 
are also presented. 

From Table I one can immediately see that the calculated masses with retardation effect 
are generally well fitted compared with those without retardation considered, and the spin 
splittings, l 3 Si — l 1 ^ are significantly improved by adding the retardation term to the 
scalar confinement potential. These conclusions obviously shed light on the usefulness of 
the Bethe-Salpeter equation in describing systems containing light quarks, however, in the 
mean time they also ask for further investigations on the interaction kernels. 

As noticed in Ref. ||, the smallness of the hyperfme splitting obtained for qQ mesons is 
due to the weakness of the binding potential in these systems. However, similarly as showed 
in Ref. [JJ this situation may also be changed after including the retardation effect in the 
interaction kernel. For demonstration, in the equal-quark-mass special case the coefficients 
for the scalar potential G s which plays the main role in setting the spin splittings in Eqs. ( [20]) 



and pip will reduce to 



n- / _ t\ 1 q ■ k — m 2 
F°s (<T; k) — — - + , (29) 

hs ^ k) ~ 2E 2 2E 2 q q 2 {E k + m y (30) 

coresponding to 0~ and 1~ mesons, respectively. It is clear that these coefficients will 
limit to m/E or it higher order while q —>■ k. On the other hand, however, the static 
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linear confining potential in momentum space, which behaves as G s (q — k) oc (q— &) -4 , is 
strongly weighted as q — * k in Eq.([20|) and Because the coefficients will diminish in the 
relativistic limit as q — > k, the strength of the confining potential would be reduced in turn. 
This is the reason which leads to the small spin splittings obtained for heavy-light-quark 
mesons. However, this depressing situation would be changed if the retardation is taken into 
account. In fact, the covariant form of confinement interaction may take the form Gs{q, k) oc 
[(q — k) 2 — (go — ko) 2 ]~ 2 and in the on-shell approximation, q 2 = m 2 + q 2 , k 2 , = m 2 + k 2 , it 
becomes 

G s (q,k) oc (~2m 2 q + qk-q-k)- 2 (31) 

in the high energy limit, i.e., p,k 3> m. We can see that with the retardation effect the 
scalar interaction Gs{q,k) is heavily weighted as q and k are co-linear (q \\ k), whereas 
the static linear potential only peaks at q = k. This indicates that the former is weighted 
in a much wider kinematic region than the latter, especially for systems including light 
components. As a consequence, the wave functions in coordinate space would tend to be 
short ranged and hence the magnitudes of the wave function at the origin, ip(0), would 
increased. Because to leading order in 1/m 2 expansion the hyperfine splitting is proportional 
to |^(0) | 2 , the Hermitan square of the wave function of meson at the origin, the splittings 
would enlarge as well. The above analysis indicates that in the equal-quark-mass situation 
the modified effective scalar interaction will not be weakened too much as q and k approach 
to be parallel, and this is just due to the retardation effect. In our opinion, the difficulty 
in the reduced Salpeter equation with the static scalar confinement is probably due to the 



improper treatment that the confining interaction is purely instantaneous JT| fl2|| . 

In practice, for the constituent quark model, which is essentially used in the present work, 
the equal mass and the on-shell approximation maybe not a good simplification for the qQ 
systems, but in any case the analysis given above is qualitatively correct, which is supported 
by the numerical results listed in Table I. And, it is obvious that the new procedure gives 
a much better fit to data than the previous one, in especially the spin splittings. The 
results obtained about the mass spectrum, particularly the s d system, seems still not fully 
satisfactory, but it is noted that our results listed are just a schematic ones. A fine tunning 
may lead an improvement on them. 

In conclusion, we have extented our previous study in clarifying the problems pointed 
out by Durand et al. for the static scalar confinement in reduced Salpeter equation to 
the heavy-light-quark systems in this paper. The conclusion remains the same that the 
"intrinsic flaw" of the Salpeter equation with static scalar confinement could be remedied 
to some extent by taking the retardation effect of the confinement into consideration. A 
fit for mass spectrum of S-wave heavy-light-quark systems is obtained by using the scalar 
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linear confinement potential with retardation and the one-gluon exchange potential. Result 
shows that a great improvement is achieved on fitting the data with retardation taken into 
account. Although the on-shell approximation may not be a rigorous treatment here, the 
qualitative feature of the retardation effect is still manifest. Nevertheless, it is still premature 
to assess whether or not the quark confinement is really represented by the scalar exchange 
of the form of (f? — £>o) -2 , as suggested by some authors and being used here, as the dressed 
gluon propagator to implement quark confinement. Therefore, further investigations on this 
subject are necessary. 
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Table I 



Calculated mass spectrum of qQ states using reduced Salpeter equation with retardation for scalar 
confinement. The Experiment data are taken from Ref. [0 



state 


Quarks 


uata ^iviev ) 


r It 1 ^ivie V ) 


ri/rror ^ivie v ) 


r It 11 ^ivie V ) 


H/rror ^ivie v ) 


D 


oct 


OA til 


Oool 


_li no 


OZOo 


-Zl 


D + 


cs 


1969 


2097 


+ 128 


1946 


-23 


Dt + 


cs 


2112 


2148 


+35 


2094 


-19 


D + 


cd 


1869 


1983 


+114 


1862 


-7 


D* + 


cd 


2010 


2010* 





2003 


-7 


K° 


sd 


498 


743 


+ 245 


652 


+154 


K*° 


sd 


892 


870 


-22 


898 


+6 


D*+ - Df 




144 


51 


-93 


148 


+4 


D* + - D + 




141 


27 


-114 


141 





K*° - K° 




394 


127 


-267 


244 


-150 



"Results without retardation. Results with retardation. *Used to fix the parameters. 
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